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TURBULENT TRANSPORT MEASUREMENTS WITH A 


LASER DOPPLER VELOCIMETER 

Robert V. Edwards , JohnC. Angus , 
and John W. Dunning, Jr. 

Lewis Research Center 
ABSTRACT 

The power spectrum of phototube current from a Laser Doppler 
Velocimeter operating in the heterodyne mode has been computed. 

The spectrum is obtained in terms of the space -time correlation func- 
tion, G(Ar, t), of the fluid. The spectral width and shape predicted by 
the theory are in agreement with experiment. For normal operating 
parameters the time -average spectrum contains information only for 
times shorter than the Lagrangian -integral time scale of the turbu- 
lence. To examine the long-time behavior, one must use either ex- 
tremely small scattering angles, much -longer -wavelength radiation 
or a different mode of signal analysis, e.g. , F. M. detection. 

1. Introduction 

In a Laser Doppler Velocimeter (LDV) laser light is scattered 
from particles suspended in the fluid. The power spectrum of the 
scattered light contains information on the motion of the particles 
and hence information on molecular and turbulent transport processes. 
A schematic view of the measurement region is shown in figure 1. 

The light received at the detector is a linear sum of electric 
fields scattered from each particle in the illuminated region. The sig- 
nal at the detector can be written 
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f(t) = A cos [w Q t + (p(t)] (1) 

where A is an amplitude, is the laser frequency, and <p( t) is the re- 
sultant phase. 

<p(t) is the phase that results from properly adding the phase received 
from each of the particles. This phase is determined by the positions of 
the particles and the scattering angle, 6. For the situation shown in fig- 

■fVi 

ure 1, the phase of light scattered from the i Ln particle is 



iHx, + c 

X 1 


(2) 


■fVi 

where is the X coordinate of the r particle, X the laser wave- 
length and a is a constant determined by the optics: at = 2 sin (e/2). 

The resultant phase, <p, and amplitude, E are obtained by summing the 
contributions, E i exp (i0.) from the individual illuminated particles. This 
is most easily seen graphically by a plot in the complex plane, figure 2. 

If the relative particle positions are fixed and if the suspending fluid 
is convecting at some constant velocity v in the X direction, the signal 
f(t) is 


• f(t) = A(t) cos 




t + C 


(3) 


The term i2ira/X) v is the frequency change usually interpreted as the 
Doppler shift. It is indeed proportional to the velocity and if this were 
all that happened, life would be simple indeed! However, as we men- 
tioned, only the particles in the illuminated region contribute to the sig- 
nal. If the particles in the illuminated region are removed and another 
set substituted by convection, the phase of the new signal is different, 
since the new particles do not occupy the same positions as the old set. 
Further, since the suspended particles are randomly placed in any real 
system, the new phase cannot be predicted solely with knowledge of the 
old phase; i.e. , the two phase terms are uncorrelated. Therefore, in 
a real system with a finite illuminated volume, the signal (eq. (3)) only 
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stays correlated with itself at time t + t as long as some particles present 
in it at t are still present at time t + t. If the particles maintain their rel- 
ative positions, this time is on the order of ct/v where a is the character- 
istic dimension of the illuminated region in the flow direction. 

The power density spectrum of a random signal is the Fourier transform 
of the autocorrelation of the signal. In this case, the autocorrelation is zero 
for all times greater than a/v so that the resulting spectrum has a width on 
the order v/a . To a first (but good) approximation this spectrum is cen- 
tered at 2nav/\, so that a measurement of the spectrum can yield v and or. 

Until this point, the particles have been assumed not to change relative 
position and that the flow was steady. If we now allow the particles to move 
relative to each other (diffuse), another phenomenon appears. In undergoing 
molecular or turbulent motion, the particles lose correlation with them- 
selves, which means that the signal loses correlation with itself. This con- 
tributes to the spectrum of the received signal. The interpretation of this 
contribution of the spectrum in terms of the stochastic properties of tur - 
bulent fluid motion is the purpose of the remainder of this paper. 

2. Theory 

The instantaneous A.C. heterodyne current from the photodetector for 
a Laser Doppler Velocimeter may be written [1], assuming a constant ref- 
erence beam 

r 

r V' iK?[r (t)l 

i(t) °c E g (t)E r (t) = vconst. R e ' p [ r n (t)] (4) 

n 

where oc implies proportionality and where r (t) is the time dependent 
position of the n L particle in the system measured with respect to a 
laboratory frame. The vector K is the scattering vector and is fixed by 
the optics. P[ r n (t)] is the weighting function for the E fields seen by the 
photodetector. The weighting function, which can be complex, is deter- 
mined by the optics and can be computed from first principles. It is es- 
sentially the product of the E field amplitude weighting function for the 
two beams that define the sample volume. The term ’’const. ” contains 
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terms describing the scattering efficiency of the particles, the quantum 
efficiency of the photodetector, etc. , which are unimportant for the pres- 
ent discussion. 

By the Weiner -Khinchin theorem, the power density spectrum of the 
photodetector current qan be written 


I(K,co) 



(5) 


where 


R^(K, = const. i(t)i(t + t) (6) 

is the autocorrelation of the current. The overbar denotes a time average. 
This now may be written in terms of the particle positions, 



Stationarity has been assumed in eqs ? (7) and (8). In all but fluid-particle sys- 
tems of very high number density, the relative particle positions can be 
taken as random. Therefore the second term in equation (8) is vanishingly 
small compared to the first term. Note here, that because of the average, 
the second term tends to zero even for the limiting case of very few par- 
ticles in the sample volume at one time. 
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Equation (8) can now be rewritten, 


= const. Re^~ 
n 


-iK- Ar (t) 

e P[r n (o)] P[r n (o) + Ar n ( T )] 


where 


( 9 ) 


Ar n (r) = r n (7) - r n (o) 

This expression can be formally evaluated in terms of two probability den- 
sity functions: 

(1) that there is a particle within the volume element r, r + dr, 


£ dr (10) 

N 


(2) that a particle in the volume element r, r + dr moves a distance 
Ar, in time r, 


G(Ar, r;r)d Ar (11) 

The autocorrelation of the photodetector current is, from equations (9), 
(10), and (11), 


R- = const. Re p 



-iK- Ar 


G(Ar, r;r)P(r)P[r + Ar(T)]dr d Ar 


( 12 ) 


3. Laminar Flow 

Consider a system in steady, straight laminar flow. In this case, 
for the scattering particles normally used in LDV systems, 

G(Ar, r, r) = 6[Ar - v(r)r] (13) 


i.e. , the particles track the flow. Equation (12) becomes 

,-iK- v(r)T. 


R^ = const. Re p / e' 


P(r)P[r 4- v(r)r]dr 


(14) 
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The dependence of v on r takes account of gradients in the velocity. 
Detailed calculations of the spectrum for laminar pipe flow have been 
presented in a previous paper [1]. A typical spectrum is shown in fig- 
ure 3. This spectrum demonstrates the effects of both the finite illu- 
iminated region and of a gradient in the velocity. Here v Q is the veloc- 
ity in the center of the sample volume and z Q /R is the fractional dis- 
tance from the center of the pipe . 

4. Turbulent Flow 

We define the average velocity v by 

rs * J 

v(r) = / — G(Ar, r;r)d Ar (15) 


For convenience, Ar will be defined in terms of the mean velocity v, 
and Ar the deviation from the mean motion 

Ar = Ar + vt (16) 


The computations for this section will be done for v * v(r). This is 
not a necessary condition, but it aids in the clarity of the presentation. 
Equation (12) now becomes 


R.. = const. 

li 



G(Ar + v T)e 


-iK- (Ar+vr) 


P(r)P(r + Ar +vr)dr d(Ar + vr) 

(17) 


In view of equation (15) we define 


G(Ar, r)d Ar = G(Ar, r)d Ar 


(18) 


where G is the probability that a particle moves Ar, Ar + d Ar with 
respect to the mean motion vr in time r. This definition becomes clear 
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if one notes that the probability that a particle moves Ar in time t is 
identical to the probability that it moves Ar with respect to the mean 
flow in time r. 

Equation (17) becomes 


R-. = const. Re j G(Ar, r)e _iK ‘ Ar e _iK ’ VT Q 0 (Ar + vr)d Ar (19) 


where 


Q 0 = / P(r)P(r + Ar + vr)dr 


(20) 


Equation (19) can always be rewritten using the Fourier convolution and 
shift theorems as 


R^ = const. 



r)e 


-iK'- vr 


Q o (K - K’)dK* 


( 21 ) 


where Q q is the Fourier space transform of Q q and} G is the Fourier 
space transform of G. 

R.. can be shown to tend to zero with increasing r for two reasons 

( 1 ) 


lim 

7“-*-oo 



-iK’ - vr 


f(K')dK' 


= 0 


( 22 ) 


for any square integrable f(K’), 

(2) G(K’, t) can be shown to be a decreasing function of r for any 
K’, in particular K, so that 


lim 


G(K’,t)Q 0 (K' - K)dK T 


= 0 


(23) 


For the first mechanism, we have demonstrated in a previous paper [1], 
that a characteristic decay time is 
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o 

where a is characteristic dimension of P in the mean flow direction. 
For the second mechanism, the characteristic decay time can be taken to 
be when 


G(K, r)/G(K, 0) e" 1 / 2 (24) 

i.e., when the standard deviation [aK 2 (t)3 ^ 2 , of G is approximately 
equal to K. But from the theory of Fourier transforms we know that 


[aK 2 (t 3 1/2 [Ar 2 (7)] 1/2 = 0(1) (25) 

2 1/2 

where (Ar ) ' is the standard deviation of G(Ar, r). 

The characteristic decay time, t for the second mechanism 
therefore occurs when 


K[Ar 2 (T Q )] 1/2 = 0(1) (26) 


In other words, it is the time when the average scattering particle has 

diffused a distance 1 /Ik| with respect to the mean motion. Since |k| 

for a typical LDV is on the order of 1 0 cm , the second mechanism 

21/2-4 

causes significant decay of when (Ar ) ' & 10 cm. This length 

is very small compared to the normally encountered Lagrangian integral 
length scales, which are on the order of the shear dimension of the 
flow system. Therefore, no matter what value Tq takes, R~ decays in 
at most the time it takes the average particle to wander a distance with 
respect to the mean flow that is very small compared to the Lagrangian 
integral length scale A^. This means that the integrand involved in the 
evaluation of R.. is significantly different from zero only for times that 
are short compared to the Lagrangian integral time scale, 

In this case, G(Ar, r) can be approximated by an asymptotic expres- 
sion for small t/t^. 


5(Ar - v’r)F(v T ;r)dv’, 



G(Ar, t) 


(27) 
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where F(v’) is the probability density function for a turbulent fluctuation 
v’. 

In this section, several simplifying assumptions are made: 

(1) G(Ar, t) is isotropic 

(2) the flow is one dimensional 

(3) the vector K lies in the mean flow direction 

These assumptions are made only for clarity of presentation and do not 
represent a restriction on the method of calculating the spectra. 

If the G(Ar, t) from equation (27) is substituted into equation (19), 
and that used in equation (5), we find after some rearrangement 

I(K,a>) = const. ^ — 1— ■ F(v?)0. f Kv - - + . ^ V ' W (28) 

^Jv’+v ^ V V' + V / 

where v* = oo/K’ - v. The integration variable, v’, in equation (28) can 
now be interpreted as the fluctuating velocity. This becomes clear when 
one recognizes that the term 

is the spectrum generated by a flow with a constant velocity v + v’ (see 
ref. 1). Equation (28) gives the spectrum as a weighted sum of spectra 
generated by each velocity present in the sample region. Individual par- 
ticles can change their velocities during their traverse of the sample vol- 
ume; but since f Q « they do not change their velocity significantly 
over, t q , the decay time for R^. 

For P = 1, a very large sample volume, 

I(K w) = const. I F ( w " K ' v ) (30) 

K \ K / 


K(v + v T ) - 


V + V’ 


(29) 


In this case the spectrum has the shape of the probability density function 
for the fluctuating velocity with a scale factor of K. 



If K is sufficiently small, G(K, t) and hence for turbulent flow, 
do not decay until times long compared to r^. In this case the observed 
particle motion takes on a Brownian motion character and G(Ar, r) can be 
approximated by (4) 

G(Ar, t) » 1 e~ Ar2 / 4eT (31) 

(4776 t) 3/2 

where e is the eddy diffusivity. Using equations (31), (21), and (5), the 
spectrum of the LDV becomes 

I(K, to) = const. / ¥—1 — Q»(K’-K)dK’ (32) 

J (K ,2 e) 2 + (u> - K’ • v) 2 

For equation (32) to be valid K must be less than l/A^. Since is 
typically greater than 1cm., it appears infeasible to measure eddy dif- 
fusivities by the above technique. However, the theory presented here 
also applies to ultrasonic scattering where small K vectors are possible. 
A check of the small K limit predicted here is feasible by ultrasonic 
techniques. For the Brownian motion of small particles, where the core- 
relation length is much smaller than 1/K, the results have been confirmed 
[2,3]. 


5. Calculation for Real Systems 

Consider a three dimensional system with a sample volume having a 

Gaussian intensity with RMS widths a , a , a in the three orthogonal 

x y Ci 

directions and with orthogonal turbulence intensities, 
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The moments of the spectrum can be calculated from the three dimensional 
form of equation (28). When the mean flow is in the X-direction, and when 
K is in the X-direction, the mean and second central moment of the spec- 
trum are 


M x = Kv 


VAR = Kv’ 2 + —+ — 


t 2 t ,»2 „,2 


x 


4a„ 4a' 


4°y *°l 



(33a) 


(33b) 


(33c) 


These formulae are independent of the particular form of the fluctuating 
velocity probability density F(v’). The v /4a terms have been discussed 
in reference [1] and are the finite sample volume contribution to the width of 
the spectrum. The other terms are the broadening due to the turbulence. 
One will notice that the variance is not simply the sum of the square of the 
finite volume broadening and the turbulence broadening. It contains, in ad- 
dition, turbulent finite volume broadening (e.g. , v? /4a 2 ) due to the RMS 
velocity fluctuations carrying the scattering centers across the sample vol- 
ume. Typically, a: /a < a /a « 1 and (v’Vv 2 )*/ 2 <0.05, so the ef- 

x z x y 

feet of y and z fluctuations is to increase the finite transit time term by 
less than 3 percent. Since, in a well -designed experiment, the finite transit 
time broadening is small compared to the turbulence broadening, the effect 
of y and z fluctuations is negligible. Therefore, a good approximation to 
the variance of the spectrum is 


VAR = K 2 v^ 2 



(34) 
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If the K -vector is in the y direction 


M 1 = 0 


( 35 ) 


VAR = K 2 v’ 2 + — 

y 9 

4 °x 


(36) 


The LDV is therefore capable of measuring each component of the fluctuating 
intensity independently. 

Near the centerline of turbulent channel flow, the gradient in the mean 
velocity is small. Further F(v') can be approximated by 


F(v’)dv’ = 


1 


,,. 2 - 


exp - 


'V' ~ V * 2 ’ V* 2 

x y z 
■ + J - ■ + — — 


27 r vj 2 vj 2 v: 2 


X y z 
Using equations (37) and (28) 


2v^ 2 2v’ 2 2v’ 2 

x y z 


d y x dVy dv^, (37) 


l(K,u>) 


const. 





< 2 K 2 

e x / x exp- 


v + v 

X X 




idv’ 

X 


(38) 


Since |v j » jv* f one has 



r . 2 


V x 


A 

exp 

2v’ 2 


X 



2or x (Kv i + Kv x " 

i . 


dv^ (38b) 
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l(K, do) 


const. 


/ 

r -2 ■ i 


/ 

\ 2jr 

Ii + 2KW 


2 j - 

V 

2 - x 




1 

!_ \ 


(co - Kv x ) 2 



(39) 


The heterodyne turbulence spect rum in this case is approximated by a 
Gaussian with variance (ir/4 a + K v’ ). Figure 4 is a turbulence spectrum 
taken at the channel center. The solid line is a Gaussian fitted to the data. 

Near the center of a channel, the dependence of the turbulence intensity 
on position is not strong, so that over the sample volume, v’ , can be con- 
sidered constant and the gradient in the mean velocity can be handled as in 
reference [1], Near walls where large gradients in t he mean velocity may 
be present and a strong positional dependence of v' 2 may be also present, 
a more complex model for G(Ar, t) must be used. This will be treated in 
a later paper . 


6 . Conclusion 

The power spectrum of phototube current for a Laser Doppler Velo- 
cimeter (LDV) operating in the heterodyne mode has been computed for the 
case of turbulent flows. It is shown that for normal operating parameters 
the spectrum contains information only on the short-time behavior of the 
fluid motion. To examine the long-time behavior, i.e. , times greater than 
the Lagrangian -integral time scale of the turbulence, one must use ex- 
tremely small scattering angles. 

It should be noted that this paper describes the power spectrum of a 
signal averaged over long times. This includes situations in which there 
may be a very few or very many particles in the sample volume at one 
time. We have not considered the case of F. M. detection. With F. M. 
detection it is possible to extend the time and length scale of the LDV 
measurement beyond those indicated in this paper. 
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Figure 1. - Sample region. 
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Figure 3. - Spectra observed with laminar pipe flow showing 
effects of large sample volume in radial (z) direction. In 
each case the sample volume was centered at 0. 1 cm from 
centerline of 1 cm pipe, the characteristic dimension, o x , 
in the flow direction was approximately 0.0045 cm and the 
velocity at the tube center was 10 cm/ sec. Note that K- v 0 
is identical in each case. 



FREQUENCY, kHz 

Figure 4. - Spectrum obtained from turbulent 
flow of water. Parameters were; v « 288 
cm/ sec, Re - 85 700, 9= 19.8°, o x = 10.4 pm. 
Calculated percent turbulence was 2.95 per- 
cent. 
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